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We consider chains of one-dimensional, piecewise linear, chaotic maps with uniform slope. We study the diffusive behaviour of 
an initially nonuniform distribution of points as a function of the slope of the map by solving the Frobenius- Perron equation. 
For Markov partition values of the slope, we relate the diffusion coefficient to eigenvalues of the topological transition matrix. 
The diffusion coefficient obtained shows a fractal structure as a function of the slope of the map. This result may be typical 
for a wide class of maps, such as two-dimensional sawtooth maps. 
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The study of simple models for non-equilibrium pro- 
cesses in statistical physics has been one of the central 
themes in the theory of chaotic dynamical systems [|l|,||,H . 
A great deal of work has been done to describe the large- 
scale motion in systems of independent particles each 
moving under the action of relatively simple maps, oper- 
ating at discrete intervals of time. For one-dimensional 
sinusoidal, or piecewise difFerentiable maps, a variety of 
diffusive-like and ballistic-like behaviour has been stud- 
ied y,|[^|J3]- For two-dimensional, conservative Hamil- 
tonian maps, parameter dependent momentum-diffusion 
coefficients have been computed, often by a combination 
of numerical and analytical methods which explore the 
phase space structure of the dynamical system p,H,nO[. 
Recently, Gaspard and coworkers have established an 
explicit connection between fundamental quantities of 
dynamical systems, such as the Kolmogorov-Sinai en- 
tropy and Lyapunov exponents, and transport coeffi- 
cients |0,0,0 14 1 . Related connections between trans- 



port coefficients and Lyapunov exponents have been dis- 
cussed for non-equilibrium systems with thermostats [Ea . 
There is also a close connection of the work described here 
to that based on periodic orbit expansions for transport 
coefficients [|7|,p^. 

In this letter, chains of piecewise linear, one- 
dimensional, chaotic maps 



Xt + 1 = [Xt] + rUaiXr) = Ma{Xr) 



(1) 



will be considered, where r is the discrete time variable, 
[xr] the largest integer smaller than x-r, ma(xr -I- 1) = 
TT^aixr) represents a periodic function, and a stands for 
the control parameter, which is the slope of the map. 
We consider a chain of maps ma{xr) with chainlength L. 
The absolute value of the slope is assumed to be uniform 
and its logarithm is equal to the Lyapunov exponent. 
We assume the maps are expanding, i.e. that \a\ > 1. 



Following the approach in ||llj,|l^,|13|,y^ , we describe a 
method by which the diffusion coefficient for this class of 
maps can be computed for a broad range of parameter 
values. The method will be illustrated by the map 
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as sketched in Fig. H, which has been introduced and 
discussed in [p|,E|Jl^]. We find that the diffusion coeffi- 
cient for this map shows a very rich fractal structure as 
a function of the slope. 

To describe the dynamical behaviour of an arbitrary 
initial density for a set of particles on some interval of the 
line — cx) < x < cx), we will need the Frobenius- Perron 
equation, supplemented by boundary conditions. The 
Frobenius-Perron equation is given by 



(3) 



Pr+i{x) = / dy priv) 5{x - Ma{v)) 



where Pt{x) is the probability density for points on 
the line, and Maijj) is the map under consideration. 
We suppose that the motion takes place on an interval 
< a; < L, and we impose periodic boundary conditions, 
i.e. Pt(0) — Pt{L) for all r, or absorbing boundary con- 
ditions Pt{x) = for X = 0, L for all r |17]. Next we use 
the argument of Gaspard and coworkers | Il| , |l2| , [l3|jl4| ] to 
relate the eigenmodes of the Frobenius-Perron equation 
to the solution of the diffusion equation 



dn{x,t) d'^n{x,t) 



dt 



dx^ 



(4) 



where n(x, t) is the macroscopic density of particles at 
a point X at time t, and £) is a diffusion coefficient. If 



for large L, and large r, the first few eigenmodes of the 
Frobenius-Perron equation are identical to those of the 
diffusion equation, the diffusion coefhcient can be ob- 
tained by matching eigenmodes in an appropriate scaling 
limit. More explicitly, for periodic boundary conditions 
n(0,i) = n{L,t), one can easily see that for large times 
n(x,t) has the form 



all of length 1/2, as illustrated in Fig. |l|. Then M is a 
2L X 2L matrix of the form 



n{x, t) = const. + Aexp {-D{At:'^ / L^)t ± i[2'K / L)x) 



(5) 



Consequently, if one can find a solution of Eq. (H), for 
large L and r, in the form of 

p{x, t) = const. + A' exp (— 7p(a)T ± i{27r/L)x) , (6) 



one can relate the decay rate 7p(a) to D by 
D{a) = lim {L/2T:f jp{a) . 

L— >oo 



(7) 



For absorbing boundary conditions one relates the dif- 
fusion coefhcient to the escape rate from the system 
by an equation similar to Eq.(n). The escape-rate for- 
malism of chaotic dynamics shows that the escape rate 
from a system with absorbing boundaries is equal to 
the Lyapunov exponent minus the Kolmogorov-Sinai en- 
tropy for particles trapped within the system whose tra- 
jectories lie on a fractal repeller mU. For the case of 
maps with slopes of uniform magnitude considered here, 
the Komolgorov-Sinai entropy for the fractal repeller 
hKsio.) is identical to the topological entropy of points 
on the repeller and it can be computed from 7p(a) as 
hKsia) = \oga-jp/A + OiL-^) |l|. 

The use of maps with uniform slope is not an essential 
ingredient in the calculation of D{a) described below, 
which can be applied to more general linear maps. The 
main idea is that the Frobenius-Perron equation can be 
written as a matrix equation whenever the parameters 
of the map are such that one can construct a Markov 
partition of the interval (0, L), which has the property 
that partition points get mapped onto other partition 
points by the map Ma{x) [Q. In a related context, these 
partitions have been discussed in ||l^ . For such values of 
a, Eq.(H) can be written as 



Pr+l = (l/|a|)^^Pr 



(8) 



where p^ is a column vector of the probability densities 
in each of the Markov partition regions at time r, and 
M is a topological transition matrix whose elements Mij 
are unity if points in region j can be mapped into region 
i, and are zero otherwise. 

As a simple example we consider the form of the matrix 
M when a — 3, the map Ma{x) is given by Eqs.(|l|),(||), 
and periodic boundary conditions are used on an interval 
of length L. In this case the regions of the partition are 
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(9) 



\0 1 ••• 1 1 / 

In the limit t — > cxd, for any L, and any "Markov par- 
tition" value of a, the Frobenius-Perron equation can 
be solved in terms of the eigenmodes of M for any ini- 
tial value po {^) which is uniform in each of the Markov 
partition regions. For periodic boundary conditions, 
M is always a (block)circulant [gOJ, the largest eigen- 
value of M is precisely \a\, and the corresponding eigen- 
mode is a constant, representing the equilibrium state. 
The rate of decay to equilibrium, 7p(a), is obtained as 
7p(a) = log(|a /xi), where xi is the next largest eigen- 
value of A/ [n3|. Analytical expressions for D{a) can be 
derived for all integer values of a > 2. For even integers, 
the results of Grossmann and Fujisaka M are recovered, 
D{a) ~ (l/24)(a — l)(a — 2), and for odd integers we 
find D{a) = (1/24) (a^ - l). To obtain D{a) for a gen- 
eral Markov partition value of a, one can use computer 
methods §^. 

Fig. g (a) shows the results for the diffusion coefficient 
of the dynamical system Eqs.(|l|),(||) for values of a in the 
range 2 < a < 8. In Fig. g (b)-(d), we present magnifica- 
tions of three small regions in this interval [2^] . One can 
see clearly that D{a) has a complicated fractal structure 
with regions exhibiting self-similarity. In Fig. @, we show 
an enlargement of the region for 2 < a < 3. The dashed 
line is the prediciton of D{a) for a simple random-walk 
model suggested by Schell, Eraser and Kapral Q. Note 
that the model correctly accounts for the behaviour of 
D(a) near a = 2. The wiggles in this graph can be un- 
derstood by considering the transport of particles from 
one unit interval to another. These regions are coupled 
to each other by turnstiles, where points in one unit in- 
terval get mapped outside that particular interval into 
another unit interval. As in the case of two-dimensional 
twist maps, such as the sawtooth map, these turnstiles 
are crucial for large-scale transport MV^ . 

The region 2 < a < 3 can be analyzed by studying the 
interaction of turnstiles |24|. One can recognize three 
distinct series of values of a, each of which provides a 
cascade of apparently self-similar regions of decreasing 
size, as the limits a ^ 2 or a ^ 3 are approached. To 
understand these series, consider the trajectory of a point 
that starts just to the left at cc = 1/2. The first iterate 
of X = 1/2 is in the second interval, (1, 2). The series a 
values of a are defined by the condition that the second 
iterate of a; = 1/2 is at the leftmost point of the upward 



turnstile in the second interval (1, 2) (a = 2.732), or that 
the third iterate is at the corresponding point in the third 
interval (a = 2.920), etc. The numbers on the graph re- 
fer to the number of intervals the image of a; = 1/2 has 
travelled before it gets to the appropriate point on the 
turnstiles. Series f3 points are defined in a similar way, 
but they are allowed to have two or more internal reflec- 
tions within an interval before reaching the left edge of 
a turnstile. Series 7 points are defined by the condition 
that some image of a; = 1/2 has reached the rightmost 
edge of an upward turnstile (i.e. some point x = n + 1/2, 
where n is an integer), and consequently an increase in a 
will lead to a decrease in D{a). These cascades provide a 
basis for a physical understanding of the features of D{a) 
in this region: Particles leave a particular unit interval 
through a turnstile and undergo a number of iterations 
before they are within another turnstile. Whether they 
continue to move in the same or the reverse direction 
at the next and later turnstiles is a sensitive function of 
the slope of the map. Thus the fractal structure of the 
D{a) curve is due to the effects of long-range correlations 
among turnstiles and these correlations lead to changes of 
D{a) on an infinitely fine scale. A similar argument can 
be employed to explain, at least qualitatively, the fractal 
structure oi D{a) for higher values of the slope, although 
more work needs to be done before a full understanding 
of this curve is obtained p3] . 

We conclude with a few remarks: (A) Our results ap- 
pear to be the first example of a system whose diffu- 
sion coefficient has an unambiguously fractal structure. 
We suspect that similar results obtain for all other one- 
dimensional, piecewise linear maps pq|, which might be 



of interest, e.g., for chaotic scattering l27l, as well as for 
transport in maps of more than one dimension, such as 
sawtooth maps. We note that oscillations of the diffusion 
coefficient with respect to an appropriate control param- 
eter, which could be a field strength, have already been 
found in the standard [g| and the sawtooth map [Q . (B) 
It is not known whether this fractal structure persists for 
smooth maps where the function Ma(x) is C^, or where 
the map contains some randomness. (C) We have numer- 
ical evidence that the Markov points are dense for a > 2, 
and we believe that our results give the full structure of 
the D{a) function. Nevertheless, it would be valuable to 
have a mathematical proof. 
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FIG. 1. Illustration of the dynamical system Eqs. ([lj),(0) 
for a particular slope, a — 3. The Markov partition given 
by the dashed grid leads to the construction of the transition 
matrix in Eq. (Mj. 



FIG. 2. Diffusion coefficient D{a) computed for the dynam- 
ical system Eqs. ([l|),(p|) and some enlargements. Graph (a) 
consists of 7908 single data points. In graph (b)-(d), the dots 
are connected with lines. The number of data points is 476 
for (b), 1674 for (c), and 530 for (d). 



FIG. 3. Enlargement of the region of slope a < 3 with 
the solution for a simple random walk model (dashed line) 
and labels for the values which are significant for "turnstile 
dynamics" (see text). For some points, the turnstile coupling 
is shown by pairs of boxes. The graph shows 979 single data 
points. 
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